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i.  INTRODUCTION 


Recent  papers1’2’3  have  reported  the  development  and  application  of  the 
thin-layer  Parabolized  Navier-Stokes  computati onal  technique  to  predict  the 
flow  about  slender  bodies  of  revolution  at  supersonic  velocities.  Reference  3 
showed  the  technique  to  be  a  viable  computati onal  tool  for  predicting  Magnus 
effects  for  a  six-caliber  slender  shell  with  a  one-caliber,  7°  boattail  as 
verified  by  comparison  to  wind  tunnel  force  measurements.  The  results  of 
Reference  3  represent  the  first  successful  efforts  to  compute  the  Magnus 
effect  for  boattailed  shell  using  sophisticated  numerical  computati onal  tech¬ 
niques.  In  this  paper,  the  ability  of  the  PNS  computational  technique  to  pre¬ 
dict  normal  and  Magnus  forces  at  angles  of  attack  up  to  10  degrees  and  the 
effects  of  nose  and  afterbody  geometry  on  the  aerodynami cs  of  shell  are 
examined.  The  influence  of  nose  bluntness  is  determined  using  a  starting 
solution  generated  by  the  unsteady  Navier-Stokes  (time  marching)  computati onal 
technique  and  then  using  the  PNS  (space  marching)  code  to  compute  the  flow 
over  the  remainder  of  the  shell. 


II.  COMPUTATIONAL  TECHNIQUE 


A.  Space  Marching  Solution 

The  steady  thin- layer  Parabolized  Navier-Stokes  equations  can  be  written 
for  general  spatial  coordinates  £,  n,  C  as1 

A 

rvp  A  ^  A 

s  .  if  _36  _i_  as  m 

an  3c  “  at;  1  1 


where 


5 

n 

c 


£(x)  is  the  streamwise  (marching)  coordinate 

n(x,.y,z)  is  the  ci rcumferenti al  coordinate 

c(x,y,z)  is  the  normal  coordinate 


1.  Schiff ,  L.B.,  and  Steger ,  J.L.,  "Numerical  Simulation  of  Steady  Supersonic 
Viscous  Flow,"  AIAA  Journal,  Voil.  18,  No.  12,  December  1980,  pp.  1421- 
1430. 

2.  Schiff,  L.B.,  and  Sturek,  W.B.,  "Numerical  Simulation  of  Steady  Supersonic 
Flow  Over  an  Ogive-Cylinder-Boattail  Body ,  "  AIAA  Paper  No.  80-0066,  AIAA 
18th  Aerospace  Sciences  Meeting,  January  1980. 

3.  Sturek,  W.B.,  and  Schiff,  L.B.,  "Computations  of  the  Magnus  Effect  for 
Slender  Bodies  in  Supersonic  Flow,"  AIAA  Paper  No.  80-1586-CP ,  AIAA 
Atmospheric  Flight  Mechanics  Conference,  August  1980. 
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The  inviscid  flux  vectors  in  Eq.  (1)  are 


pU 

P  v 

pW  ““ 

PuU+C  p 

^  -1 
,  F  «  J  1 

puV+nxp 

A  -1 
,  G  *  J 

puW+cxP 

A  ) 

pvU 

pvV+n^p 

pvW+i^P 

pwU 

pwV+nzp 

pwW+^zp 

(e  +  p$)U 

(e  +  p)V 

(e  +  p)W 

The  vector  q(p,  pu,  pv,  pw,  e)  contains  the  dependent  variables.  Variations 
of  body  geometry  are  included  in  Eq.  (1)  through  the  presence  of  the  metric 
terms  nx> '  n  ,  nz>  etc.,  which  appear  in  the  flux  vectors  and  through  the 

Jacobian  J,  of  the  transformation  from  the  physical  space  to  the  computational 
space.  The  thin-layer  viscous  term,  valid  for  high  Reynolds  No.  flow,  is 


S  =  J 


-1 


*  (>'/3)(?xu4Hyv5H2w4k)! 

+  <“/3)(ixut+ytnzwt)sy 
»('x+ty+'z>"{  +  (W3)(txucHyvsH2wt)cz 
fUx»GyHz)  [ ( u/2 ) (u2+v2+w2 ) t  +  «Pr"1(T-l)“1(a2)t] 
+  (m/3)  (cxu  +  ?yv  t  C2“)(Cxuc  +  +  S2»c)l 


Equation  1  is  parabolic-like  with  respect  to  5  and  can  thus  be  marched  down¬ 
stream  in  the  ?  direction  from  an  initial  data  plane  (subject  to  appropriate 
body  and  free  stream  boundary  conditions)  under  those  conditions  where  the 
local  flow  is  supersonic. 

The  numerical  algorithm  used  to  march  Eq.  (1)  downstream  is  an  approxi¬ 
mately  factored,  fully  implicit,  finite-difference  scheme.  The  algorithm  is 
conservative  and  of  second-order  accuracy  in  the  marching  direction.  A  two- 
layer,  algebraic  eddy  viscosity  model4  is  included  for  the  computation  of 
turbulent  flows.  Details  of  the  Parabolized  Navi er-Stokes  assumption  and  the 
derivation  of  the  algorithm  are  included  in  Reference  1. 


4.  Baldwin,  B.S. ,  and  Lomax ,  H.,  "Thin  Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257,  1978. 
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The  computations  are  started  from  a  converged  conical  solution  near  the 
tip  of  the  projectile.  The  full  solution  is  then  obtained  by  marching  over 
the  body  in  the  streamwise  (axial)  direction.  This  marching  technique  is 
depicted  in  Figure  1  along  with  the  coordinate  system.  A  logarithmic  stretch¬ 
ing  is  used  to  achieve  adequate  grid  resolution  of  the  turbulent  viscous 
layer.  The  grid  generator  employs  an  adaptive  capability  which  insures  that 
adequate  resolution  of  the  viscous  layer  is  maintained  over  the  full  length  of 
the  model.  The  streamwise  marching  stepsize  was  adjusted  to  yield  500  to  800 
computational  steps  for  the  full  length  of  the  shell.  The  computational  grid 
consisted  of  36  stations  about  the  circumference  of  the  model  (A<j>  =  10°)  and 
50  points  between  the  body  and  the  outer  boundary. 

B.  Blunt  Nose  Solution 


The  solution  for  the  flow  over  the  blunt  nose  of  the  shell  has  been 

obtained  using  the  three  dimensional  thin-layer  Navi er-Stokes  solver  recently 
reported  by  Chaussee,  Kutler  and  Pulliam.5  This  is  a  fully  implicit  computa¬ 
tional  technique  which  solves  for  the  entire  flow  field  by  converging  to  a 
steady  flow  solution.  The  outer  boundary  for  this  code  employes  a  shock¬ 

fitting  scheme.  Since  the  solution  is  obtained  for  only  a  small  flow  field 
region  near  the  nose,  good  grid  resolution  has  been  achieved.  The  grid  used 
consisted  of  30  points  from  the  body  to  the  outer  boundary,  20  points  axially 
along  the  surface  of  the  shell,  and  36  points  circumferentially  around  the 
shel 1 . 

The  nose  region  has  initially  been  modeled  as  a  hemisphere  cone.  The 
conical  extension  is  of  sufficient  length  to  achieve  a  well  established  super¬ 
sonic  conical  flow  in  order  to  apply  the  zero  gradient  outflow  boundary  condi¬ 
tion  with  confidence.  This  modeling  of  the  nose  region  is  depicted  schematic¬ 
ally  in  Figure  2.  The  bluntness  ratio  is  defined  to  be  the  ratio  of  the 
diameter  of  the  cone  at  the  forward  tangency  point  to  the  maximum  diameter  of 
the  shell. 

In  addition  to  the  hemisphere  blunt  nose,  several  cases  have  been  run  for 

a  flattened  nose.  The  intent  is  to  model  as  closely  as  possible  the  meplate 

(flat  face)  configuration  of  typical  Army  fuze. 

The  time-dependent  thin-layer 
strong  conservation-law  form  as 


aq  +  a£ 

9t  95 


Navi er-Stokes  equations  can  be  written  in 


3n  H 


_l  dS 
Re  85 


(2) 


5.  Chaussee,  D.S.,  Kutler,  P.,  and  Pulliam,  T.H.,  "Three  Dimensional  Viscous 
Flew  Field  Program;  Part  I:  Viscous  Blunt  Body  Program  (Interim  Report)," 
A FWL- TM-81 -63- FIMG,  March  1981. 
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where 


T 


•-  t  is  the  time  (marchinq) 

5  =  S(t,  x,  y,  z)  is  the  longitudinal  coordinate 

n  =  n(t,  x,  y,  z)  is  the  circumferential  coordinate 

C  =  c(t,  x,  y,  z)  is  the  near  normal  coordinate 


AAA 

The  inviscid  flux  vectors  are  E,  F,  6  and  contain  terms  for  the  conservation 
of  mass,  momentum,  and  energy  in  the  three  coordinate  directions.  The 

vector  q(p,  pu,  pv,  pw,  e)  contains  the  dependent  variables.  The  S  matrix 
contains  the  viscous  terms  which  are  valid  for  high  Reynolds  number  flows. 

The  contravarient  velocity  components  are 


u  =  sxu  +  y  +  y 

V  =  nxu  +  HyV  +  nzw 
W  =  sxu  +  cyv  +  czw. 

Variations  of  body  geometry  are  included  in  Equation  (2)  through  the  presence 
of  metric  terms  (5  ,  n  ,  etc.)  and  the  Jacobian  J,  which  appear  in  the  flux 

vectors  and  contravariant  velocities.  A  two-layer  algebraic  eddy  viscosity 
model  is  used  for  the  turbulent  viscous  solutions. 


III.  RESULTS 

A.  Comparisons  to  Aerodynamic  Force  Measurements 

In  Reference  3  a  series  of  results  obtained  for  a  =  2°  were  compared  to 
experimental  measurements  (Reference  6)  of  aerodynamic  forces  for  2  <  M  <  4. 
The  model  configuration  for  this  study  is  shown  in  Figure  3.  The  sign  conven¬ 
tion  for  the  aerodynamic  forces  is  shown  in  Figure  4.  The  results  in 

Reference  3  indicated  that  the  thin-layer  PNS  computational  technique  achieved 
excellent  agreement  in  comparison  to  experimental  measurements  of  aerodynamic 
forces  including  the  Magnus  effect.  Comparisons  of  the  PNS  computed  results 
to  experiment  have  now  been  carried  out  for  a  <  10°.  Computed  results  were 
obtained  for  a  =  2°,  4°,  6°,  and  10°. 


6.  Nietubicz ,  C.J .,  and  Opalka,  K.,  "Supersonic  Wind  Tunnel  Measurements  of 
Static  arid  Magnus  Aerodynamic  Coefficients  for  Projectile  Shapes  with 
Tangent  and  Secant  Ogive  Noses,"  ARBRL-MR-02991,  U.S.  Army  Ballistic 
Research  Laboratory/ ARRADCOM,  Aberdeen  Proving  Ground,  MD  21005,  February 
1980  (AD  A083297 ) ." 
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Comparisons  between  computation  and  experiment  are  shown  as  a  function  of 
angle  of  attack  for  normal  force  for  M  =  2  and  M  =  3  in  Figures  5  and  6, 

respectively.  Results  for  two  boattail  confi gurati ons  are  shown.  The  nota¬ 
tion  SOCBT  depicts  the  7°,  one  caliber  boattail  model  in  Figure  3.  The  nota¬ 

tion  SOC  depicts  a  model  with  a  0°  boattail.  Excellent  agreement  is  achieved 
for  a  <;  10°  for  magnitude  as  well  as  the  nonlinear  behavior  of  the  normal 
force  coefficient. 

A  similar  comparison  is  shown  for  Magnus  force  for  M  =  2  and  M  =  3  in 

Figures  7  and  8,  respectively.  These  results  indicate  very  good  agreement 
between  computation  and  experiment  for  a  <  6° .  However,  the  computed  results 
for  a  =  10°  are  not  in  close  agreement  with  the  trend  of  the  experimental  data 
for  a  >  6°.  This  is  more  clearly  indicated  in  Figure  8.  This  discrepancy  is 
not  unexpected  since,  as  shown  in  Reference  2  where  comparisons  between  compu¬ 
tation  and  experiment  for  surface  pressure  were  made,  good  agreement  was 
achieved  at  a  =  6.3°  and  significantly  less  satisfactory  agreement  was  obtain¬ 
ed  at  a  =  10.4°.  The  discrepancy  at  a  =  10.4°  was  attributed  to  the  inability 
of  the  computational  technique  to  accurately  model  the  significant  development 
of  lee-side  vortical  flow. 

The  ability  of  the  computational  technique  to  accurately  predict  Magnus 
for  a  <  6°  does  represent  a  significant  capability.  This  demonstrates  that 
useful  engineering  results  can  be  achieved  for  highly  three-dimensional  flow 
fields  using  a  very  simple  turbulence  model. 

B.  Effect  of  Boattail  Geometry 


1 .  Scope  of  Computational  Effort. 

A  series  of  computations  have  been  accomplished  for  a  parametric  vari¬ 
ation  of  boattail  configuration.  The  geometries  are  shown  in  Figure  9  and 
include  boattail  lengths  of  one  and  two  calibers  for  boattail  angles  of  0° , 
5°,  7-1/2°,  and  10°.  This  range  of  boattail  length  and  angle  effectively 
spans  the  range  for  practical  shell  application. 

The  computations  were  accomplished  for  standard  atmospheric  and  wall 
temperature  conditions  commonly  encountered  in  projectile  firing  tests.  These 
conditions  are  summarized  in  Table  1. 
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TABLE  1.  SUMMARY  OF  BOUNDARY  CONDITIONS  FOR  PARAMETRIC  COMPUTATIONS 


MACH 

a, 

T„. 

P„» 

V 

No. 

degrees 

PD /V 

°K 

atm 

IL 

2 

2 

.19 

294 

1 

239,294,325 

3 

2 

.19 

294 

1 

239,294,325 

4 

2 

.19 

294 

1 

239,294,325 

MACH 

No.  FREE  STREAM  REYNOLDS  NO. 


2  4.53  x  102 3 * * *  7/in 

3  6.80  x  10 7/m 

4  9.06  x  107/m 


2.  Aerodynamic  Forces  versus  Axial  Position. 

The  development  of  the  normal  force  as  a  function  of  axial  position  is 
shown  in  Figure  10  for  boattail  angles  of  0°  and  7-1/2°  comparing  results  for 
one  and  two  caliber  lengths  of  the  boattail.  The  normal  force  increases 
monotonical ly  for  the  0°  boattail  in  contrast  to  the  behavior  for  the  7-1/2° 
boattails  where  the  normal  force  reaches  a  maximum  and  then  decreases  as  the 
boattail  length  increases. 

The  development  of  the  Magnus  force  as  a  function  of  axial  position  is 
shown  in  Figures  11  and  12  comparing  the  effects  of  boattail  length  (Figure 
11)  and  boattail  angle  (Figure  12).  Figure  11  shows  that  the  length  of  the 
boattail  strongly  affects  the  magnitude  of  the  Magnus  force.  Figure  12  shows 
that  the  Magnus  force  increases  monotoni cal ly  for  increasing  boattail  angle 
for  the  cases  considered  here. 

3.  Aerodynamic  Coefficients  versus  Mach  Number. 

Parametric  comparisons  for  pitch  and  yaw  plane  aerodynamic  parameters 
are  shown  in  Figures  13  and  14.  These  examples  illustrate  the  ability  of  the 
computational  technique  to  develop  data  which  reflect  the  effects  of  body  con¬ 
figuration  and  Mach  number. 

The  effect  of  boattail  angle  and  length  for  the  slope  of  the  pitching 
moment  (CG  =  3.6  calibers)  is  shown  in  Figure  13.  These  results  show  that 
both  boattail  length  and  angle  have  a  significant  effect  on  the  pitching 
moment.  The  trends  indicate  that  CM  is  increased  for  increasing  boattail 

a 

angle  and  boattail  length.  This  trend  is  accentuated  for  flow  velocities  near 
Mach  =  2. 
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The  slope  of  the  Magnus  moment  coefficient  (CG  =  3.6  calibers)  is 
shown  in  Figure  14.  The  effect  of  boattail  length  is  seen  to  strongly  affect 
the  Magnus  moment  coefficient.  The  trend  illustrated  is  that  is 
increased  as  boattail  length  and  angle  are  increased.  pa 

C.  Effects  of  Nose  Bluntness 


A  spark  shadowgraph  of  an  artillery  shell  at  a  speed  of  approximately  Mach 
=  2.25  taken  in  the  BRL  Transonic  Range  is  shown  in  Figure  15.  This  shadow¬ 
graph  shows  the  detached  bow  shock  that  occurs  for  supersonic  flow  over  artil¬ 
lery  shell.  Note  the  presence  of  a  shock  at  the  sharp  corner  of  the  flattened 
nose.  This  shock  is  caused  by  the  local  flow  separation  which  is  induced  by 
the  strong  expansion  at  the  sharp  corner. 

A  recent  wind  tunnel  test7  was  performed  in  which  aerodynamic  force 
measurements  were  made  for  spinning  models  with  sharp,  hemisphere,  and  flat¬ 
tened  nose  tips.  A  Schlieren  photograph  of  the  flow  over  these  nose  tips  at 
Mach  =  3  is  shown  in  Figure  16.  It  is  obvious  that  the  bow  shock  pattern  is 
different  for  each  of  these  configurations.  The  inner  shock  at  the  nose  tip 
that  is  very  distinct  in  Figure  15  also  occurs  at  the  sharp  corner  of  the 
flattened  nose  in  Figure  16.  This  shock  is  not  visible  for  the  hemisphere 
nose.  The  shock  waves  starting  about  one  diameter  downstream  of  the  model 
nose  are  generated  by  a  boundary  layer  trip. 

In  this  portion  of  the  study,  the  significance  of  the  technique  for  model¬ 
ing  the  nose  region  of  spinning  shell  on  the  predicted  aerodynamic  behavior  at 
supersonic  velocities  is  examined.  The  model  geometry  used  is  the  ogive- 
cyl inder-boattail  (SOCBT)  shape  shown  in  Figure  3.  Figure  2  shows  the  nose 
tip  geometries  in  detail.  Flow  field  computations  have  been  obtained  for  wind 
tunnel  test  conditions  and  for  free  flight  atmospheric  conditions.  The  wall 
temperature  boundary  condition  has  been  modeled  as  an  adiabatic  wall  for  the 
wind  tunnel  tests  and  as  a  constant  wall  temperature  equal  to  the  atmospheric 
temperature  for  the  atmospheric  flight  conditions. 

Examples  of  u-velocity  profiles  at  the  starting  plane  for  sharp  and  blunt 
noses  are  shown  in  Figure  17.  The  profile  for  the  blunt  nose  cases  are  signi¬ 
ficantly  distorted  in  the  inviscid  flow  region  compared  to  the  sharp  nose  case 
which  has  a  comparatively  flat  profile.  An  interesting  comparison  is  shown  in 
Figure  18  which  illustrates  Mach  contours  for  the  blunt  starting  solutions. 
The  more  severe  flow  expansion  around  the  flattened  nose  compared  to  the 
hemisphere  nose  is  clearly  illustrated.  Figures  19  and  20  show  the  distribu¬ 
tion  of  surface  pressure  for  windward  and  leeward  rays  for  the  two  blunt  nose 
shapes. 

Additional  examples  of  u-velocity  profiles  for  the  different  nose  config¬ 
urations  are  shown  in  Figure  21  for  a  longitudinal  station  near  the  start  of 
the  boattail.  Small  but  distinct  differences  are  apparent  upon  comparison  of 


7.  Unpublished  Wind  Tunnel  Data.  Test  performed  at  the  Naval  Surface  Weapons 
Center,  White  Oak  Laboratory,  Silver  Spring,  Maryland,  for  the  U.S.  Army 
Ballistic  Research  Laboratory. 
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these  profiles  indicating  that  the  effect  of  the  initial  conditions  does  not 
wash  out  rapidly. 

The  wind  tunnel  test7  conducted  recently  for  the  BRL  at  the  Naval  Surface 
Weapons  Center,  White  Oak  Laboratory  obtained  measurements  of  the  aerodynamic 
forces  on  the  slender  shell  with  sharp,  hemisphere,  and  flattened  nose  config¬ 
urations.  Force  balance  measurements  were  obtained  for  pitch  plane  and  Magnus 
forces.  Computational  results  for  sharp  and  hemisphere  nose  geometries  have 
been  obtained  for  comparison  to  these  data.  The  computational  and  experiment¬ 
al  results  for  Magnus  force  are  compared  in  Figure  22.  The  experimental  data 
indicate  an  increasing  Magnus  force  for  increasing  nose  flatness.  The  agree¬ 
ment  between  computation  and  experiment  for  the  sharp  and  hemisphere  nose  con¬ 
figurations  is  quite  good.  However,  there  is  substantial  disagreement  between 
computation  and  experiment  for  the  flattened  nose  geometry.  The  computation 
for  the  hemisphere  and  the  flattened  nose  yielded  virtually  identical  results. 
This  is  in  sharp  contrast  to  the  experimental  results  which  yielded  a  substan¬ 
tially  greater  Magnus  force  for  the  flat  nose  case  than  that  for  the  hemi¬ 
sphere  cap.  Further  computations  for  which  the  boundary  layer  on  the  nose  cap 
was  laminar  yielded  results  that  were  virtually  identical  to  the  results  for 
which  the  boundary  layer  on  the  nose  cap  was  turbulent.  A  possible  cause  for 
the  discrepancy  is  the  difference  between  the  flatness  of  the  nose  for  the 

computational  model  and  the  model  used  in  the  wind  tunnel  test.  As  shown  in 
Figure  2,  the  flat  part  of  the  nose  does  not  extend  to  the  full  diameter  of 
the  meplate.  Additional  computations  are  planned  in  which  the  actual  flatness 
of  the  fuze  configuration  will  be  modeled  more  closely.  It  would  also  be 

useful  to  have  experimental  data  for  a  wider  range  of  free  stream  Mach  number 
for  comparison  to  the  computations. 

An  example  is  shown  in  Fi qure  23  of  the  development  of  the  Magnus  force  as 
a  function  of  axial  position  comparing  computational  results  for  sharp,  hemi¬ 
sphere,  and  flattened  nose  configurations  for  atmospheric  flight  conditions. 
The  results  indicate  a  significant  increase  in  the  Magnus  force  for  the  blunt 
noses.  The  trend  of  increasing  Magnus  force  for  increasing  nose  flatness  is 
observed.  Examples  of  the  Magnus  force  as  a  function  of  Mach  number  for  the 

sharp  and  hemisphere  blunt  noses  are  shown  in  Figure  24.  The  blunt  case  is 

consistently  greater  in  absolute  magnitude  than  the  sharp  case. 

The  slope  of  the  Magnus  moment  coefficient  as  a  function  of  Mach  number  is 
shown  in  Figure  25  comparing  results  for  sharp  and  hemisphere  noses.  Again, 
the  blunt  nose  results  are  consistently  greater  than  those  for  the  sharp  nose. 
The  results  also  predict  a  greater  influence  of  nose  bluntness  as  the  Mach 
number  decreases.  The  behavior  of  the  Magnus  center  of  pressure  is  shown  in 
Figure  26.  These  results  indicate  that  the  Magnus  center  of  pressure  is:  (1) 
located  well  downstream  of  the  CG  location  (3.6  calibers),  (2)  weakly  sensi¬ 
tive  to  Mach  number,  and  (3)  not  highly  sensitive  to  the  nose  bluntness. 

The  slope  of  the  pitching  moment  coefficient  as  a  function  of  Mach  number 
is  shown  in  Figure  27.  These  results  indicate  that  the  pitching  moment  is  not 
sensitive  to  small  nose  bluntness  typical  of  artillery  shell. 
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IV.  SUMMARY 


A  computati onal  study  has  been  described  in  which  thin-layer  Navier-Stokes 
computati onal  techniques  have  been  employed  to  predict  the  aerodynamics  of 
slender,  spinning  shell  at  supersonic  velocities. 

Results  have  been  discussed  which  illustrate  the  ability  of  the  computa¬ 
tional  technique  to  accurately  predict  Magnus  and  normal  forces  for  angles  of 
attack  up  to  six  degrees.  Parametric  results  showing  the  effects  of  boattail 
length  and  boattail  angle  were  discussed  for  pitching  moment  and  Magnus 
moment. 

Comparisons  have  been  shown  between  results  obtained  in  which  the  nose  tip 
of  the  shell  has  been  modeled  as  a  sharp  cone,  a  hemisphere-cone,  and  a 
truncated  cone.  The  results  indicate  that  small  nose  bluntness,  typical  of 
artillery  shell,  can  significantly  increase  the  Magnus  moment.  The  results 
further  indicate  that  the  pitching  moment  is  not  strongly  affected  by  the 
small  nose  bluntness  considered  in  this  study. 
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Figure  1.  Coordinates  and  Notation. 
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Figure  2. 


Illustration 


of  Blunt  Nose  Modeling  Geometries. 
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NOTE: DIMENSIONS  ARE  IN  MILLIMETRES 

A,B  BOUNDARY  LAYER  SURVEY  STATIONS 
C  CONICAL  STARTING  SOLUTION  STATION 


Figure  3.  Model  Configuration  for  Wind  Tunnel  Test. 
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NORMAL  FORCE 


Figure  4.  Sign  Convention  for  Aerodynamic  Forces. 
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—  PNS  COMPUTATION,  SOC 

—  PNS  COMPUTATION,  SOCBT 
•  EXP,  SOC 


Figure  5.  Normal  Force  Coefficient  versus  Angle  of  Attack, 

Computation  Compared  to  Wind  Tunnel  Force  Measurements, 
M  =  2,  Re  =  2.49  x  107/m 


—  PNS  COMPUTATION,  SOC 
--  PNS  COMPUTATION,  SOCBT 

•  EXP,  SOC 

♦  EXP,  SOCBT 


Figure  6.  Normal  Force  Coefficient  versus  Angle  of  Attack, 

Computation  Compared  to  Wind  Tunnel  Force  Measurements, 
M  =  3,  Re  =  2.11  x  107/m 
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Figure  7.  Magnus  Force  Coefficient  versus  Angle  of  Attack,  Computation 
Compared  to  Wind  Tunnel  Force  Measurements,  M  =  2, 

Re  =  2.49  x  107/m,  PD/V  =  .23,  Tw  =  Taw 


—  PNS  COMPUTATION,  SOC 


Figure  8.  Magnus  Force  Coefficient  versus  Angle  of  Attack,  Computation 
Compared  to  Wind  Tunnel  Force  Measurements,  M  =  3, 

Re  =  2.11  x  107/m,  PD/V  =  .19,  Tw  =  Taw 
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Figure  9.  Model  Configurations  for  Parametric  Study 
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Figure  10.  Normal  Force  Coefficient  versus  Axial  Position,  Parametric 

Comparison,  M  =  3,  a  =  2°,  T  =  T 


Figure  11.  Magnus  Force  Coefficient  versus  Axial  Position, 

Parametric  Comparison,  M  =  3,  a  =  2°,  PD /V  =  .19, 
Tw  =  Taw,  Re  =  6.80  x  10?/m 


Figure  12.  Magnus  Force  Coefficient  versus  Axial  Position, 

Parametric  Comparison,  M  =  3,  a  =  2°,  PD /V  =  .19, 

Tw  =  Taw>  Re  =  6-80  x  lo7/m 
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Figure  13.  Slope  of  Pitching  Moment  Coefficient  versus  Mach  No.,  Parametric 
Comparison,  =  294K,  Tw  =  294K,  CG  =  3.6  calibers 


Figure  14.  Slope  of  Magnus  Moment  Coefficient  versus  Mach  No.,  Parametric 
Comparison,  =  294K,  Tw  =  294K,  CG  =  3.6  calibers 
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Figure  15.  Aerodynamics  Range  Spark  Shadowgraph  of  Artillery 
Shell  at  Supersonic  Velocity,  M  «  2.3 


26 


Figure  16.  Schlieren  Photographs  of  Wind  Tunnel  Flow  over  SOCBT  Shell  with 
Sharp,  Hemisphere,  and  Flattened  Nose  Geometries,  M  =  3, 

Re  =  9.19  x  106/m 
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Figure  17.  Lee  Side  U-Velocity  Profiles  at  Starting  Plane,  M  =  3,  a  =  2°, 
PD/V  =  0.19,  Atmospheric  Flight  Conditions,  Re  =  6.80  *  lO'/m 
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Figure  18.  Mach  Contours  of  Nose  Region  Flow  for  Hemisphere  and  Flattened 
Nose  Geometries,  M  =  3,  a  =  2°,  PD /V  =  0.19,  Atmospheric  Flight 
Conditions,  Re  =  6.80  x  107/m 
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Figure  19. 


Surface  Pressure  on  Hemisphere  Nose,  M  = 
Atmospheric  Flight  Conditions,  Re 
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Figure  20. 
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Surface  Pressure  on  Flat  Nose,  M  =  3,  a  =  2°,  PD/V 
Atmospheric  Flight  Conditions,  Re  =  6.80  x  10'/ro 
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Figure  21.  Lee  Side  U-Velocity  Profiles  at  X/D  =  4.88,  SOCBT  Model, 

M  =  3,  a  =  2°,  PD/V  =  0.19,  Atmospheric  Flight  Conditions, 
Re  =  6.80  x  107/m 
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Figure  22.  Magnus  Force  versus  Axial  Position,  SOCBT  Model,  Computation 
Compared  to  Experimental  Force  Measurements,  M  =  3,  a  =  2°, 
PD/V  =  0.19,  Wind  Tunnel  Test  Conditions,  Re  =  9.19  x  106/m 
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Figure  23. 


Magnus  Force  versus  Axial  Position,  SOCBT  Model,  M  =  3,  a  =  2°, 
PD/V  =  0.19,  Atmospheric  Flight  Conditions,  Re  =  6.80  x  107/m 
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Figure  24.  Magnus  Force  versus  Mach  Number,  SOCBT  Model,  Comparison  between 
Sharp  and  Hemisphere  Noses,  a  =  2°,  PD /V  =  0.19,  Atmospheric 
Flight  Conditions 
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Figure  25.  Slope  of  Magnus  Moment  Coefficient  versus  Mach  Number, 
Model,  Comparison  between  Sharp  and  Hemisphere  Noses, 
Atmospheric  Flight  Conditions 
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Figure  26.  Magnus  Force  Center  of  Pressure  versus  Mach  Number,  SOCBT  Model, 
Comparison  between  Sharp  and  Hemisphere  Noses,  Atmospheric 
Flight  Conditions 
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Figure  27.  Slope  of  Pitching  Moment  Coefficient  versus  Mach  Number,  SOCBT 
Model,  Comparison  between  Sharp  and  Hemisphere  Noses, 
Atmospheric  Flight  Conditions 
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LIST  OF  SYMBOLS 


a 

=  speed  of  sound 

CP 

=  center  of  pressure  for  normal  force 

CPY 

=  center  of  pressure  for  Magnus  force 

CM 

=  pitching  moment  coefficient 

CM 

pa 

=  d2Cn/ [d(-^j)  *da] ,  slope  of  Magnus  moment  coefficient  evaluated  at 
PD/V  =  0,  a  =  0 

cn 

*  Magnus  (yawing)  moment  coefficient 

CN 

=  normal  force  coefficient 

Cy 

=  Magnus  (side)  force 

D 

=  diameter  of  model 

e 

=  total  energy  per  unit  volume  of  fluid,  normalized  by  p^a^2 

A  A  A  A 

e,es,f,g 

=  flux  vectors  of  transformed  gasdynamic  equation 

J 

=  Jacobian  of  transformation  between  physical  and  computational 
coordi nates 

L 

=  reference  length 

M 

=  Mach  number 

P 

=  pressure  normalized  by  pa  2 

Pr 

=  Prandtl  number,  p  c  /< 

5  OO  p'  GO 

PD/V 

=  nondimensional  spin  rate  about  model  axis 

Re 

=  Reynolds  number,  pJJJ-/^ 

A 

s 

=  viscous  flux  vector 

u.v.w 

=  Cartesian  velocity  components  along  the  x,  y,  z  axis, 
respectively,  normalized  by  aa 

u.v.w 

=  contravari ant  velocity  components 

x.y.z 

=  physical  Cartesian  coordinate  axes 

a 

=  angle  of  attack 

Y 

=  ratio  of  specific  heats 
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LIST  OF  SYMBOLS  (Cont'd) 


K 

M 


5,n,C 


p 


=  coefficient  of  thermal  conductivity,  normalized  by  free- stream 
value  k 

00 

=  coefficient  of  viscosity,  normalized  by  free-stream  value 

*  computational  coordinates  in  the  axial,  circumferential,  and 
radial  directions  (Fig.  1) 

=  density,  normalized  by  free-stream  density  pw 


Subscripts 


aw 

w 

x 


=  free-stream  conditions 
=  adiabatic  wall 
=  body  surface  values 
=  based  on  axial  distance  from  nose 
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